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■ Abstract. 

^ . DECIGO Path Finder (DPF) is a space-borne gravitational wave (GW) detector 

(-H I with sensitivity in the frequency band 0.1-lOOHz. As a first step mission to DECIGO, 

it is aiming for launching in 2016-2017. Although its main objective is to demonstrate 

Q . technology for GW observation in space, DPF still has a chance of detecting GW 

^ ' signals and performing astrophysical observations. With an observable range up to 

C/2 . 50 kpc, its main targets are GW signals from galactic intermediate mass black hole 

^ ' (IMBH) binaries. By using inspiral-merger-ringdown phenomenological waveforms, we 



perform both pattern-averaged analysis and Monte Carlo simulations including the 
effect of detector motion to find that the masses and (effective) spins of the IMBHs 
I could be determined with errors of a few percent, should the signals be detected. Since 

GW signals from IMBH binaries with masses above lO^M© cannot be detected by 
ground-based detectors, these objects can be unique sources for DPF. If the inspiral 
signal of a 10'^ M0 IMBH binary is detected with DPF, it can give alert to the ringdown 
signal for the ground-based detectors lO^-lO'^s before coalescence. We also estimate the 
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, possible bound on the graviton Compton wavelength from a possible IMBH binary in 

uj Ccntauri. We obtain a slightly weaker constraint than the solar system experiment 
and an about 2 orders of magnitude stronger constraint than the one from binary 
_ pulsar tests. Unfortunately, the detection rate of IMBH binaries is rather small. 
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1. Introduction 

There are several ground-based interferometers aiming for the first direct detection 
of gravitational waves (GWs) from astrophysical sources such as inspirals of neutron- 
stars (NSs), black holes (BHs) and mixed binaries, and also from supernovae and 
gamma-ray bursts (GRBs). Among them, LIGO is now under upgrading phase to 
advanced LIGO [Il[2]- VIRGO is currently operating, but will soon become offline and 
will be upgraded to advanced VIRGO pi H], while GEO [5l E] will continue to take 
data. TAMA is also in upgrading phase to LCGT [71 [8]. These detectors have their 
best sensitivity around 100-lOOOHz and their lower frequency sides are limited by the 
seismic noiseq||. To overcome this limit, space-borne interferometers have been proposed. 
Among them, the Laser Interferometer Space Antenna (LISA) has been proposed with 
its optimal sensitivity around the GW frequency of ImHz [9l [10]. The expected 
targets are supermassive black hole (SMBH) binaries and white dwarf (WD) binaries. 
Unfortunately, the latter GWs may cover up other signals including a primordial GW 
background [HI [12] . As a prototype mission, LISA path finder (LPF) will be launched 
before LISA [13]. 

Another space-borne interferometer, the Deci-Hertz Interferometer Gravitational 
Wave Observatory (DECIGO) [HI [15] has been proposed in Japan. It is most sensitive 
at around 0.1-lHz. Similar interferometer, the Big Bang Observatory (BBO) has been 
suggested as a follow-on mission to LISA [T6l WJ\ [18] . Thanks to the high frequency 
cutoffs of WD /WD binary signals at / ~ 0.2Hz [19], the main target of DECIGO 
and BBO is the primordial GW background (PGWB). There is a possibility of NS/NS 
foreground signals masking PGWB but it has been shown that BBO is expected to 
have enough sensitivity to subtract sufficient amount [13 [20]. DECIGO may have to 
improve its sensitivity by 2-2.5 times in order to achieve this goal [21]. They have other 
interesting scientific objectives including direct detection of the accelerating expansion of 
the universe [H] , performing the precision cosmology [18] , revealing the thermal history 
of the early universe [22l [23] and the formation process of SMBH [24], verifying the 
primordial BH as the dark matter candidate [251 126] , and probing alternative theories 
of gravity [27] and the size of extra dimension [28] . 

DECIGO Path Finder (DPF) is planned as a first step mission for DECIGO, 
hopefully launched in 2016-2017 [291 [30]. Its main goals are to test the key technologies 
for the space mission and to perform observations of GWs and Earth gravity. As for 
the first space-borne GW detector, the torsion-bar type space antenna called SWIM has 
already been launched in 2009 and successfully performed its observation run [31]. It 
has completed its mission and is currently under the phase of data analysis. Recently, 
the improved version of the torsion-bar antenna (TOBA) has been proposed by Ando et 
al. [31] . In this paper, we focus on the scientific significances of GW observations with 
DPF. 

The main source for DPF is the intermediate mass BH (IMBH) binaries, where an 
I Initial LIGO had a cutoff frequency at 40Hz. 
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IMBH refers to a BH having a mass of (See e.g. Refs. [351 [33] for the studies 

of detecting GWs from IMBH binaries with LISA.) There is plenty of evidence that there 
exist stellar-mass BHs and SMBHs, but there is no direct detection of individual IMBHs 
and their existence is still disputed. One of the possible evidence for their existence is the 
discovery of the ultraluminous X-ray sources (ULXs) [3l] whose luminosity exceeds the 
Eddington luminosity of a stellar-mass BH, though their dynamical friction timescales 
are longer than the ones for SMBHs. Other possible evidence includes the radial mass- 
to-light ratio [35] and the surface brightness profiles [36l |37] of globular clusters, stellar 
proper motions near their centers [381 ESj and observations of several millisecond pulsars 
in the galactic globular cluster [IQ]. IMBH detection is important since it will give us 
clues for the formation mechanism of SMBHs [H]. Since the observable range of DPF 
for these sources are within 50kpc, DPF targets are IMBH binaries in our Galaxy. There 
are possibilities that IMBHs exist at the center of globular clusters and massive young 
clusters (see Refs. [Ml W2\ for reviews on IMBHs) and they may form binaries [131 HI]. 

We perform Fisher analysis using the spin-aligned phenomenological inspiral- 
merger-ringdown waveform [l5] to estimate how accurately we can measure binary 
parameters if the signal has been detected. We perform both pattern-averaged analysis 
and Monte Carlo simulations, taking the motion of DPF into account for the latter. We 
also comment on the possible joint search with DPF and ground-based detectors. DPF 
may be able to give alert to the latter about 10 mins before the coalescence. Also if 
only the ringdown signal is obtained with the ground-based interferometers, DPF data 
including inspiral and merger information may help in distinguishing the signal from 
noises. Furthermore, by combining the data of DPF and the ground-based detectors, 
the ringdown efficiency e^d [IHl |47] can be determined, which cannot be measured from 
the ringdown signal alone due to the degeneracy against the distance to the source. 

Following Ref. [H], we also consider the possible constraint on the graviton 
Compton wavelength \g from the DPF observations of IMBH binaries in our galaxy 
(see Refs. [^ 150] for the reviews on massive gravity theories). We compare our result 
with the current constraints on A^, (i) \g > 2.8 x 10^''cm [Slj from the solar system 
experiment and (ii) > 1.6 x lO^^cm [52] from binary pulsar observations. DPF 
constraint would be important since it is the bound in the strong-field regime, whereas 
the current constraints mentioned above have been obtained both in the weak-field 
regime. 

Unfortunately, the expected detection rate of IMBH binaries with DPF is rather 
low. We found that it is 10~^yr~^ for comparable-mass IMBH binaries and 10"^- 
lO^^yr"^ for intermediate-mass ratio inspirals (IMRIs) using advanced DPF. 

DPF has an ability to detect the gravitational wave background (GWB) with the 
energy density of flow ^ However, a stronger bound has already been set from the 

Big Bang Nucleosynthesis (BBN). For the observation of the gravity of the Earth, DPF 
is expected to perform complementary operation compared to other missions such as 
CHAMP, GRACE and GOCE [53]. (DPF can still observe GWs since geogravity noise 
dominates other noises only below / = 0.03Hz.) However, we do not consider these 
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issues further in this paper and stick to the observation of GWs from IMBH binaries. 

This paper is organized as follows. In Sec. |2l we review current observational results 
of IMBHs and explain the formation mechanisms of IMBHs and their binaries. In Sec. El 
we first review the design and concepts of DPF. Then, we introduce its noise sensitivity 
and compare it with the ones of the ground-based detectors. We derive beam-pattern 
functions for 1-armed interferometer and explain how to take the effect of DPF motion 
into account. In Sec. HI we explain our numerical setups and show the results of both 
pattern-averaged and Monte Carlo simulations. We point out possible joint searches 
with the ground-based detectors in Sec. [5] and show the possible constraints on the 
graviton mass with DPF in Sec. [6l In Sec. [3, we calculate the expected detection rate 
of DPF and summarize our work in Sec. [HI We take c = G = 1, Hq = 72km/s/Mpc, 
Vim = 0.3 and f^A = 0.7. Throughout this paper, we neglect the eccentricities of the 
binaries. 

2. Intermediate-Mass Black Holes 

The main target of DPF is IMBH binaries in our Galaxy. In this section, we briefly 
review the current observational implications for the existence of IMBHs and the 
formation mechanisms of IMBHs ^\\ [Ml [M] and IMBH binaries [13]. 

2.1. IMBH Observations 

It has been discovered that there exist compact objects with masses larger than ?)Mq 
from the careful measurements of the radial velocity of the companion stars. These 
objects are considered as the stellar-mass BHs which formed from the core collapses 
of massive stars. Also, it is likely that there exist SMBHs at the centers of galaxies. 
However the formation process of SMBHs has not been established and one possible 
way is the collisions of IMBHs with masses IO^-IO^Mq. 

One type of evidence for the existence of IMBHs is the existence of ULXs whose 
luminosity exceeds lO^^erg/s, the Eddington luminosity of a IOM0 BH (see e.g. Ref. [M] 
for a review on ULXs observations). This implies that ULXs have masses greater than 
IOMq. On the other hand, they exist several hundred pc away from the galactic centers 
on average [55]. Since the dynamical friction timescale of SMBHs is smaller than the 
Hubble time, these BHs are considered to have sunk to the galactic centers by now. 
This indicates that the masses of ULXs are smaller than lO^M©. These facts lead to 
the conclusion that they are likely to be IMBHs. (There are alternative explanations 
for ULXs other than IMBHs such as standard stellar-mass BHs with jets or relativistic 
beamings [56].) 

Another type of evidence is obtained from the radial mass-to-light ratio M /L which 
is estimated from the observed profile of the line-of-sight velocity dispersion a. Gerssen 
et al. [35] showed that there may exist an IMBH at the center of the galactic globular 
cluster M15. Noyola et al. [36] observed the surface brightness profile of the globular 
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Table 1. The distances and possible IMBH masses at the centers of galactic globular 
clusters. 



NGC 


distance 


(total) BH mass 


No. 


(kpc) 


(Mo) 


5139 (w Cen.) 


4.8 [eg 


(3.0-4.75) X 10^ [Sg 






<1.2xl0^ [381139] 






(1.3-2.3) xlO^ [37] 


6388 


10.0 [63] 


5.7x10^ [60] 


6715 (M54) 


26.8 [63] 


9.4x10^ [HI] 


6752 


4.0 [63] 


2.0x10^ [10] 


7078 (M15) 


10.3 [63] 


3.2x10^ [35] 



cluster oj Centauri, one of the largest and most massive galactic globular clusters. They 
claimed clear rise in the velocity dispersion a towards the center and excluded the 
"no BH" case with more than 99% confidence level. However, Anderson and van der 
Marel [38l [39] analyzed a catalog of 10^ stellar proper motions near the center of uj 
Centauri but did not find any significant rise in the density profile towards the center. 
They obtained an upper bound on the mass of IMBH at the center if it exists. Recently, 
Miocchi [37] analyzed the surface brightness profile of u Centauri and claimed that the 
mass of dark object at the center should lie in the range 1.3 x IO^Mq < M < 2.3 x IO^Mq. 
Also, significant core rotations have been observed in the clusters mentioned above 
and 47 Tuc [571 EH] which may be signs for the existence of IMBH binaries [59] . The 
possible masses of IMBHs at the centers of the galactic clusters and their distances are 
summarized in Table [H 

Also, observations of several millisecond pulsars in the galactic globular cluster 
NGC 6752 indicate that there exists a (l-2)xlO^M0 object within the inner 0.08pc of 
the cluster, assuming that the positive spin derivatives of the pulsars are due to the 
acceleration by the cluster gravitational potential well [ID]. Furthermore, Gebhardt et 
al. [6l] and Gerssen et al. [35] found that the rotational speed of the center of M15 is 
comparable to the velocity dispersion but A^-body simulation predicts no rotation in the 
cluster core when there is no massive compact object in the cluster. Interestingly, this 
rotation can be explained when there exists a (20 + 1O^)M0 BH binary at the center of 
the cluster with its orbital separation a ~ 10~^pc [34] • IMBHs are interesting for cluster 
dynamical evolutions and as the sources for GWs. 

2.2. IMBH Formation 

There are several IMBH formation mechanisms proposed in globular clusters. The 
first one is the repeated hardening of relatively massive BH and lOM© stellar-mass BH 
binaries via three-body interactions and their mergers due to gravitational radiation. 
Since the escape velocity of the typical globular cluster is about 50km/s, massive BHs 
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should have masses larger than 50Mq in order to stay in the cluster [HS] • If initial stars 
obey the Salpeter mass function, about 10~^ of stars have masses larger than SOM©, 
yielding tens of massive BHs in the cluster. At the center of the cluster, there are 
comparable amounts of stellar-mass BH and main sequence stars with their number 
densities n = lO^pc"'^. Within a Hubble time, initially a 50Mq BH increases its mass 
to about IO^Mq. 

The second mechanism is large BHs directly capturing stellar-mass BHs. If a stellar- 
mass BH passes near a larger BH sufficiently close, it gets bounded and merges quickly 
via gravitational radiation. There is also a formation mechanism proposed by Miller 
and Hamilton [66], which uses Kozai resonance. A binary-binary interaction results in 
a hierarchical triple system. When certain conditions are realized, the Kozai resonance 
works, making the eccentricity and the inclination of the inner binary to oscillate. If 
the maximum eccentricity becomes around unity, inner binary ends in a quick merger 
due to the gravitational radiation. 

In young clusters with the ages less than a few times 10'' yr, the most massive 
stars are still on the main sequence. These will sink to the center of a cluster 
and a core collapse takes place only among the massive stars. The resulting high 
central density leads to the runaway collisions of stars and a very massive star (VMS) 
formation [671 [68| |69 | [70]. This will eventually collapse to IMBH. This runaway growth 
occurs generically in clusters with collapsing times shorter than 3 Myr [71]. IMBHs may 
also be produced from the gravitational collapse of the population III stars [34j . 

2. 3. IMBH Binary Formation 

Binary stars affect the dynamical evolution of globular clusters via dynamical encounters 
and binary star evolutions. Ivanova et al. [72] performed full Monte Carlo N-body 
simulations and showed that in order to explain the current binary fraction, the initial 
fraction needs to be almost 100%. Then, Giirkan et al. [13] performed Monte Carlo 
simulations of runaway stellar collisions in young dense clusters, and found that when the 
initial fraction is larger than 10%, two VMSs larger than IO^Mq form as consequences 
of the runaway collisions. First, binary-binary induced runaway collisions take place 
at off-center of the cluster, producing the first VMS. Since binary-binary interactions 
destroy binaries, soon the core binary population is depleted, leading the cluster cores 
to collapse. During this phase, the second VMS forms via runaway collisions induced 
by single-single scatterings. These VMSs are expected to collapse to IMBHs, eventually 
forming an IMBH binary. 

Another possibility for IMBH binary formation is via cluster mergers. A large 
number of young binary clusters have been observed in the Magellanic Clouds [73]. If 
they contain IMBHs at their centers, they are expected to form binaries after their 
mergers. 



GW observations of galactic IMBH binaries with DPF 



7 



3. DECIGO Path Finder 

3.1. Design and Concepts 

DPF is a prototype mission of DECIGO [29|, [30] to test the advanced key technologies 
of DECIGO such as (i) a precise position measuring system with Fabry-Perot (FP) 
cavity, (ii) a highly stabilized laser source and (iii) the drag-free control system which 
shields external forces caused by solar radiation and residual atmosphere. The weight 
of the satellite is about 350kg and it will be orbiting the Earth at the Keplerian velocity 
with an altitude of 500km for 1 yr observation. It contains 2 mirrors forming a 1-arm 
interferometer with armlength 30cm. This interferometer includes a FP cavity with a 
finesse of about 100 and the laser is emitted from a highly stabilized source with an 
output power of lOOmW at a wavelength of 1030nm. Currently, the FP cavity has not 
been tested in space and DPF is expected to have better sensitivity than LPF which 
uses a Mach-Zender interferometer [13]. DPF provides new possibilities for a precise 
position measurement and high-stabilized laser in a space environment. 

Each mirror is placed inside a module called housing. There are electrostatic-type 
local sensors on the frame of the housing, which are used for measuring the relative 
positions of the mirrors and the frame. The common motion signals of 2 mirrors are 
fed back to the satellite position using thrusters (drag-free control) while the differential 
motion signals are sent to the actuators on the frame of the housing to stabilize the 
FP cavity. The drag-free controls have been performed by several satellites such as 
TRIAD-I and Gravity Probe-B satellites. LPF will operate them at the Lagrange 1 
(LI) point where the gravitational environment is stable. On the other hand, DPF will 
demonstrate it in an Earth orbit. This will open a new window for future space missions. 

The main GW sources for DPF are IMBH binaries in our Galaxy. DPF observation 
is important because observational frequency around 0.1-lHz cannot be reached by 
ground-based interferometers and LPF. Also the development of data analysis technique 
for DPF has a significant meaning since these data are expected to be more complicated 
than the one from ground-based interferometers due to the satellite orbital motion and 
the effects of the Earth. DPF can also measure the gravity of the Earth with comparable 
sensitivity to other space missions currently operating and can provide complementary 
observation to others [53] . 

3.2. Noise Spectrum 

When we assume that the noise n{t) is stationary, the noise spectral density S'„(/) can 
be defined as 

{n*{f)n{n)=Kf-n\s^{f\ (1) 

where n{f) denotes the noise in Fourier domain and (. . .) represents the expectation 
value. The root noise spectral density ^/SJJ) of DPF is depicted in Fig. [1] as the (red) 
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Figure 1. The root noise spectral density Sn{f) of DPF is shown as the (red) thick 
solid curve, together with the one for adv. DPF (high frequency part is limited by the 
shot noise rather than the laser frequency noise) as the (green) thin solid curve. We 
also show the adv. LIGO low frequency cutoff lOHz as the (magenta) thick dotted line 
and the one for LCGT down to IHz as (blue) thin dotted curve. The sky-averaged 
amplitude of a IO^Mq equal-mass binary in w Centauri is depicted as (light blue) thick 
dotted-dash curve and the one of a equal-mass binary in NGC 6752 is drawn as 

(black) thin dotted-dashed curve. Here, we assumed that the observations are carried 
out from / = 0.03Hz to mergers. We also assumed that these BHs have dimensionless 
effective spin angular momentum x — 0.2. 




thick solid curve. It is given as [74] 

«/)^i^(^)" + 4.0xlO-»H.-. (2) 

where the first term corresponds to the acceleration noises due to the collision of residual 
gas molecules, the Earth gravity, the magnetic fields of the satellite and the thermal 
radiation of the housing, while the second term represents the laser frequency noise. 
We set the lower frequency cutoff at / = 0.03Hz due to the Earth gravity. It may be 
possible to improve the sensitivity (we call this the "adv. DPF") whose noise spectral 
density is given as [75] 

This time, the high frequency part of the sensitivity is limited by the shot noise. Its 
root noise spectral density is shown as the (green) thin solid curve in Fig. [1] The 
(magenta) thick dotted vertical line at / = lOHz corresponds to the adv. LIGO lower 
cutoff frequency while the (blue) thin dotted vertical one corresponds to the root noise 




Figure 2. We introduce the detector coordinate {x,y,z} with its origin placed at 
the center of the DPF arm. The x-axis coincides with the detector arm while z-axis 
is orthogonal to this x-axis and exists on N-x plane, with TV denoting a unit vector 
pointing towards the source. The angle of N measured from x is denoted as 6 while 
the polarization angle ijj is defined as the angle from one of the polarization axes p to 
n which is a unit vector made by projecting z onto the plane perpendicular to TV. h 
is orthogonal to TV and lies on x-z plane. 



spectral density of LCGT with / = IHz cutofl§|. The latter can be expressed as [77] 

/ f \ -40 / p \ -14 

\LCGT,low _ A A in-25 J \ i i c i n-32 / / 



^n(/)^^^'''°" = 4.4xlO-^^(^j +1.8xl0-^^( ^) . (4) 



3. 3. Beam-Pattern Functions and the Effect of Detector Motion 

In this subsection, we first derive the beam-pattern functions for the 1-armed 
interferometer. Let us consider the situation shown in Fig. [2l The detector arm lies on 
the X-axis and N denotes the unit vector from the detector to the source. The z-axis 
lies on the x-N plane and is orthogonal to the x-axis. L is the unit orbital angular 
momentum vector of the binary and n is a unit vector constructed by projecting z onto 
the plane perpendicular to N . We define p = which forms one of the polarization 

principal axes. The other principal axis q is defined as q = — • is the angle of N 
measured from x and the polarization angle ip is the angle of n measured from p. 
By using the polarization tensors defined as 

Hab = PaPb - qaQb , H^f, = PaQb + QaPb , (5) 

the beam-pattern functions for a two-armed interferometer with arm directions x°' and 
y°' under x^-jja = can be defined as [THl [79] 

^ - ry") , ^ l^.ix'^x' - ry") ■ (6) 

§ It is likely that adv. LIGO and adv. VIRGO also have non-zero sensitivity at / < lOHz, but for this 
frequency rage, LCGT has better sensitivity compared to them. For comparison, Einstein Telescope 
(ET) is planned to have a sensitivity of y/tjj) = 5 x lO^^ijj^-i/a at / = IHz [^. 
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The beam-pattern functions for a one-armed interferometer can be obtained by setting 
y« = (0, 0, 0) as 

F+{d, i)) = ^H+x^x'' = ^ sin^ 6 cos2^, (7) 

{d, i)) = ^H^f,x''x'' = ^ sin^ 9 sin 2ip. (8) 

Their sky-averaged values of them are {F^'^) = (F^^) = 2/15. In Fig.[Tl we also show the 
sky-averaged GW amplitudes of a 10"* Mq equal-mass IMBH binary in u Centauri and a 
IO^Mq equal-mass one in NGC 6752 as the (light blue) thick dotted-dash and the (black) 
thin dotted-dashed curves, respectively. We use the phenomenological inspiral-merger- 
ringdown hybrid waveforms (explained in Sec. 14.21) for a circular spin-aligned binary to 
estimate these amplitudes. We assume the dimensionless effective spin parameter x of 
these BHs as x = 0.2, where x is defined a^ 
1 + 6 1-6 

X = Xi + X2 • (9) 

Here, 6 = (mi — m2)/M and Xi = Si/mf with and 5*^ being the mass and the spin 
angular momentum of z-th BH, respectively, x ranges from -1 to 1, and for the equal 
spin case, x becomes X = Xi = X2- X = 0.2 is selected as the most probable from Fig. 1 
of Ref. [80] (see also Ref. [81] for a recent measurement of the spins of the BHs at the 
centers of globular clusters). It can be seen that the GW frequency of IMBH binary 
in u Centauri is too low for the ground-based interferometers. Therefore GW signals 
of this kind become unique sources for DPF. (Higher harmonic signals of ringdown [82] 
may be detected with LCGT if it is sensitive enough down to IHz.) On the other hand, 
the one in NGC 6752 can be detected with both DPF and the ground-based ones. It 
may be possible to perform joint searches between these detectors (see Sec. [5] for more 
details) . 

In Fig. [31 we show the (sky-averaged) observable range of DPF with the signal-to- 
noise ratio (SNR) threshold set to pt^r = ^ as the (red) thick solid curve (see Eq. (fT6ll 
for the definition of SNR). The (green) thin solid curve corresponds to the one with 
adv. DPF. We also show the one of adv. LIGO and LCGT with the same curve as 
in Fig. [TJ The (black) dashed horizontal line at = 8.5kpc represents the galactic 
center. We also plot the possible GW sources in the globular clusters shown in Table [H 
assuming that they consist of equal-mass IMBH binaries. 

Next we introduce the celestial coordinate {x, y, z} shown in Fig. HI x-axis points 
the vernal equinox and z-axis points the north celestial pole and is orthogonal to the 

II In general, each BH spin has three degrees of freedom. However, since we here assume that each BH 
spin is (anti-)paraUcl to the orbital angular momentum, there is only one degree of freedom left (i.e. 
its magnitude). Furthermore, the dominant contribution of the spin in the waveform appears with the 
combination shown in Eq. ([9]) [45]. The hybrid waveform that we discuss in Sec. 14.21 can be accurately 
described with only one (effective) spin parameter x [45] . 

% This choice of pthr is rather optimistic and the false alarm rate may not be so small, but wc think 
that pthr = 5 is the minimum value that we can use to claim the detection (or at least its possibility) 
of GW signals [75]. 
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Figure 3. This figure shows the observable range for each detector with SNR thresliold 
Pthr = 5. The (red) thick solid curve represents the pattern-averaged one for DPF while 
the (green) thin solid curve shows the one with adv. DPF. Other 2 curves correspond 
to the same detectors as in Fig. [TJ We also plot the masses and distances for assumed 
equal- mass IMBH binaries in the galactic globular clusters shown in Table [TJ 





ecliptic 



vernal 
equinox 



Figure 4. We introduce the celestial coordinate {x, y, z} with x-axis pointing the 
vernal equinox and the z-axis is orthogonal to the celestial plane and is pointing the 
celestial North Pole. DPF orbits the Earth with its plane always facing towards the 
Sun and its arm pointing the Earth. It has an inclination angle id = 97.3° against the 
celestial plane. The magnified picture of the DPF orbital plane is shown at the top 
left panel. The angular velocity of the Earth we is measured from the x-axis while the 
one of the DPF ujd is measured from the axis obtained by projecting the z-axis onto 
the DPF orbital plane. 



celestial plane. We express {9,ip) in terms of the source direction {9s,4>s) and the 
direction of the orbital angular momentum (^l, 4>l) below. 

First, 6 is given as cos 6 = x-N with N = (sin cos 0s, sin sin 0g, cos ^s)- (Here, 
we assumed that distances to sources are much larger than lAU.) As shown in Fig. El 
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DPF orbits the Earth anti-clockwise seen from the Sun, with its arm pointing towards 
the Earth. The orbit is solar synchronous and it has orbital inclination of zd = 97.3°. 
At the time t, x is approximately given as 

^ — cos (y9E(t) cos COS^D — sin </?e(^) sin ^ 

x= — sin(/?E(t) cos<y5D(t) cos^D + cos(y9E(t) siny9D(^) , (10) 

COS (/9d (i) sin ZD y 

with v'e(^) = i^E^ + V^Eo and ^puit) = + V^do- Here, Ue = 27r/lyr and Wd = 
\/ (Re+500km)3 each represents the angular velocity of the Earth orbiting the Sun and the 
detector orbiting the Earth, respectively, and (pEo and ip^o denote the position of the 
detector at t = 0. When deriving Eq. ( ITOl) . for simplicity, we assumed that the Earth 
orbits in the celestial plane rather than in the ecliptic plane. The difference is negligible 
as long as the source is situated sufficiently far away from the Sun {D^ ^ lAU). Then, 
cos 9 becomes 

cos 9 = cos ipD (t) sin i^ cos 

- [cos{v2E(i) - 0s}cosv2d(^) coszd + sin{(^E(^) - 0s} sin(y9D(i)] sin^s. (11) 

Next, the polarization angle ip is estimated as cosip = h ■ p. By using L = 
(sin^L COS0L, sin^L sin (^L, cos ^l) and n = (cot 6')iV— (sin6')"^a;, it becomes 

. N X L 

cos w = n ■ — — 

\N X L\ 

cos^s{- cos{v2e(^) - 0l} sin<^D(^) + sin{v9E(^) - 0l} cosv9D(t) cosio} sini9L 
+ {cos{(pE{t) - 0s} sin (puit) cos 

— cos(/)D(^){sin{<y5E(t) — 0s} cos^L cos^D + sin{0L — 0s} sin^L sinin}} sin^s 

1/2 



'1 — COS 6'){1 — (sin^s sin6'L cos(0L — 0s) + cos6's cos6'l) } 



(12) 



4. Binary Parameter Estimation 

In this section, we perform Fisher analysis to estimate how accurately we can determine 
the binary parameters if the signals are detected. 



4.I. Fisher Analysis 

If a signal has been detected with a SNR much greater than 1 
accuracy of the parameter 9^ is given as 

A^' = x/(F^, 

where the Fisher matrix Fjj is defined as 



the determination 



(13) 



dh 



dh 



(14) 
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with the inner product 



(A|i?)^4Re/ df^^^£^. (15) 







SNR is defined by using this inner product as 

p ^ (16) 

4-2. Waveform Modeling for IMBH Binaries 

In this paper, we use the phenomenological inspiral-merger-ringdown waveform 
developed by Ajith et al. ^51 [86] . They basically matched the Taylor Tl post- 
Newtonian (PN) inspiral waveform [87] with the merger and ringdown waveform 
obtained via numerical simulations. Then, they parameterized this hybrid waveform 
with binary parameters. The Fourier transform of the waveform can be expressed as 
h{f) = A(/)e**''-^^ where the amplitude A{f) and the phase \&(/) are givenj^n Eq. (1' 
of Ref. [15]. For the coefficient C that appears in the amplitude, we use 



, 5 M^/^ 

where Apoi(t(/)) is defined as 

ApoiW)) = V(l + cos2 tyF+{e, V)2 + 4(cos2)2Fx {e, 7/;)2 (18) 

with i representing the inclination angle of the binary. The sky-averaged values of the 
beam-pattern functions are {F'^'^) = (F^^) = 2/15 which yield the sky-averaged value 
of C as (C) = 3^^2/3^^ • 

When performing Monte Carlo simulations, we use the waveform 

h{f) = ^(/)e~'['^(-/')+'^p°i(*(-^))+'^D(*(-^))], (19) 

where </3poi(^(/)) is the polarization phase given as 

, ,,,,, (l + cos^i)F+(t(f)) 
cos[^poi(t(/))] = ^ fiff^^ ^ (20) 

2cosiFX(t(/)) 

and ifj){t{f)) is the Doppler phase defined as 

^dW)) = 27r//2sin^scos[0(t(/)) -0s] (22) 
with R = lAU. 

+ For the time t{f), we take up to 3.5PN order. This can be obtained by integrating the inverse of 
df /dt given in Arun et al. |89| with respect to / and adding the spin-orbit and spin-spin couphngs that 
appear at 1.5PN and 2PN order, respectively (see e.g. Ref. [88]). 
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4-3. Numerical Setups 

In this subsection, we explain how we performed the Fisher analysis numerically. We 
set the cutoff frequencies of DPF as (/low, /high) = (0.03Hz, lOOHz). Then, we take the 
integration range of (/in, /hn) with 

/in = max{/iow, /lyr}, /fin = min{/high, /s}. (23) 
/lyr is the frequency at 1 yr before coalescence, which is given as [HSl EH] 



/lyr = 7.38 X 10" 



Hz. (24) 



Here Ai = Mrf^^ is the chirp mass where M is the total mass and 77 = /i/M is the 
symmetric mass ratio with the reduced mass /i. /s is the ringdown cutoff frequency given 
in Ref. [15]. We impose a prior distribution on x as |x| < 1- Following Berti et al. [88] . 
we perform the numerical integration with the Gauss- Legendre routine GAULEG |90j . 
This quadrature uses the zero points of the n-th Legendre polynomials as the abscissas 
and the integrand can be calculated exactly up to (2r;,-l)-th order. For taking the 
inversion of the Fisher matrix, we use the Gauss- Jordan elimination [90]. In order 
to make sure that the inversion has been performed correctly, we first normalize the 
diagonal components of the Fisher matrix to unity. Then we perform the inversion and 
convert it to the inverse of the original Fisher matrix (see Appendix C in Ref. |91j). 
We multiply the original Fisher matrix with the numerically obtained inverse matrix 
and see how close the result is to the identity matrix 5ij to check our inversion scheme. 
Among the analyses below, the Monte Carlo simulations for the massive gravity theories 
(Sec. 16.31) is expected to be the most difficult in taking the inverse of the Fisher matrix 
since it involves largest number of parameters (11 in total). Even in this case, we found 
that the difference in each component between the numerically calculated and the exact 
identity matrices is of (9(10~^) at most. (10~^ is the criterion used in Berti et al. [88].) 



4.4- Pattern- Averaged Analysis 

In this subsection, we show the results for the pattern-averaged (sky-averaged) analysis. 
Here, we have 6 binary parameters in total as 

61 = (lnM,lnr/,x,to,0o,/^L) (25) 

with to and 00 each represents the coalescence time and phase, respectively, and 
Dl denotes the luminosity distance. We estimate Fisher matrix by using the 
phenomenological waveform explained in Sec. I4.2l and calculate the determination errors 
of the binary parameters with DPF. We assumed that IMBH binaries at -Dz, = 5kpc 
contain equal-mass BHs with x = 0.2. We also set to = 0o = 0. In the panel (a) of 
Fig. [5l we show the estimation errors of InM (red thin solid). In 77 (green thick solid), 
X (blue thin dotted) and ItiDl (magenta thick dotted) against the total mass M. The 
corresponding SNRs are shown in the panel (b) as red thick solid curve. For an equal- 
mass IMBH binary in u Centauri (dotted-dashed vertical line in the panel (a)), we see 
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Figure 5. (a) DPF Parameter determination accuracies of InM (red thin solid), In 77 
(green thick sohd), x (blue thin dotted) and InD^ (magenta thick dotted) against 
total mass M for equal-mass binaries of x = 0.2 at Dl = 5kpc, obtained from pattern- 
averaged analysis. The vertical dottcd-dashed line represents a possible IMBH binary 
in oj Centauri. (b) The SNRs for equal- mass binaries ofx = 0.2atDL = 5kpc against 
total mass M with DPF (red thick solid) and with adv. DPF (green thin solid). The 
horizontal dashed line corresponds to the detection threshold pthr — 5. 



that parameters have only a several % errors. Therefore DPF may accurately determine 
the binary parameters if the signals are detected. If we use adv. DPF, roughly speaking, 
the parameter determination accuracies scale linearly with SNRs shown as a (green) thin 
solid curve in the panel (b). For example, the ones for a (10^ -|- 10^)Mq IMBH binary 
in u Centauri improves roughly by a factor 7. The horizontal dashed line at p = 5 
corresponds to the detection threshold. Since SNRs are not so high, these are only 
approximate estimations [92] . 



4-. 5. Monte Carlo Simulations 

Next, following Refs. [HH [911 [27], we performed Monte Carlo simulations. This time, 
we have 10 parameters as 

61 = (lnM,lnr/,x,to,0o,/^L,^s,0s,^L,0L). (26) 

We consider lO^M© equal-mass IMBH binaries with x = 0.2 at 5 kpc. For the direction 
of the sources (^s, (^s); we take the one of uj Centauri which has galactic latitude b and 
longitude / as (bj) = (14.97°, 109.10°) [63|. This can be converted into (^s,0s) by the 
relations explained in Appendix A Then we randomly generate 10^ sets of cos i9l in the 
range [-1,1] and 0l, '^eo and ipEo in the range [0,27r]. For each set, we estimate Fisher 
matrix and obtain the probability distribution of the determination errors at the end. 

In Fig. [6l we show the probability distributions for the determination errors of the 
binary parameters using DPF. In the panel (a), we show the ones for InM, In 77 and x- 
We see that these parameters can be determined with a several % errors. In the panel 
(b), we show the ones for In Dl and the angular resolution AQs with the latter defined 
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Figure 6. (a) Histograms showing the parameter determination accuracies of InM 
(red thin sohd), In 77 (green thick soUd) and x (blue dotted) using DPF obtained from 
Monte Carlo simulations. We assumed 10^A/q equal-mass IMBH binaries in a; Centauri 
with X = 0.2 and random orientations for the orbital angular momenta, (b) Histograms 
showing the parameter determination accuracies of In (magenta dotted-dashed) and 
fis (blue sohd) using DPF. 



as 

Afi, ^ 2vr| s\nh\^jT~^^~~^T^^. (27) 

Here, Sj^ corresponds to the inverse of the Fisher matrix. It can be seen that these 
parameters are poorly determined. The angular resolution can be roughly estimated as 

Aa~(-5j .2.4^ (j^^^j [-n^) ' 

where A is the wavelength of GWs and D is the effective size of the detector which 
we take as the diameter of Earth. This estimate shows a good agreement with 
our Monte Carlo result. In Dl is determined from the GW amplitude but since the 
angular resolution is 0(1), AlnD^ also becomes 0(1). Therefore, if we take source 
directions and orientations into parameters, it is difficult to determine (^s, 0s, ^l, 0l) 
and InD^. However, thanks to the weak degeneracies between these parameters and 
(InM, ln?7, x), the latter set can still be determined fairly accurately. Since ground- 
based interferometers are not sensitive for GW signals from IMBH binaries with BH 
masses IO^Mq, DPF will perform unique operations for these sources. 

5. Jointed Searches with Ground-Based Interferometers 

From Fig. [H we see that the signals from IO^Mq equal-mass IMBH binary in NGC 
6752 can be detected with both DPF and the advanced generations of ground-based 
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10^ 



Figure 7. Accumulated sky-averaged SNRs with DPF (red thick sohd) and with 
adv. DPF (green thin sohd) against the time to coalesce for a 10'^ Af© equal-mass 
IMBH binary with x = 0.2 in NGC 6752. We also show the one for a (10-^2 x 10^)Mq 
binary in NGC 6752 with adv. DPF (blue dotted). The horizontal dotted-dashed line 
corresponds to the detection threshold of pthr — 5. 

interferometers. The (red) thick sohd curve in Fig. [7] shows the accumulated SNR of 
this signal with DPF against the time to coalesce. DPF starts detecting this signal at 
lO^s before coalescence and when t 4 x lO^s, the accumulated SNR goes beyond 5. 
When t = lis, the GW frequency reaches / = IHz and LCGT may start detecting the 
signal. This means that DPF and adv. DPF may be able to claim detection about 7 
mins before coalescence if data analysis can be performed simultaneouslj0- In order 
to accomplish this goal, DPF needs to send the data to Earth at least every 7 mins. 
Then DPF can give alert to the ground-based detectors so as to make sure that they 
are operating and get ready for the detections. By combining the data of DPF and 
the ground-based detectors, we can see the whole history of late inspiral, merger and 
ringdown of the binary. (See a related work by Amaro-Seoane and Santamaria [93] on 
possible joint observations of IMBH binaries with LISA, ET and adv. LIGO.) Also, DPF 
data may help in confirming the actual GW signal detected by the ground-based ones 
when only the ringdown phase has been detected since they are difficult to distinguish 
from the noise^. If we use adv. DPF, the time when SNR reaches 5 becomes t ~ lO^s. 

* The number of GW cycles for this signal would be roughly 0.03Hz x lO^s ^ 3 x 10^. As for the 
data analysis, since the duration of the signal is comparable to (or slightly longer than) the orbital 
period of the satellite around the Earth (100 mins), we need to calibrate the Doppler-shifted signal 
and take the change in the position of the satellite into account. However, these can be performed 
beforehand so that the fundamental part of the data analysis should be almost the same as the usual 
matched filtering analysis. It would be something similar in between the chirp signal analysis on the 
ground-based detectors and continuous wave search from pulsars |75| . 

jj If the galactic IMBH ringdown signals are detected with the ground-based interferometers, the SNRs 
would be considerably large. Therefore, these signals may be confirmed only by taking cross-correlations 
of the ground-based detectors. 
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In order to give an earlier alert, it is important to reduce the acceleration noise rather 
than the shot noise. For example, we found that if we improve the sensitivity of DPF 
by a factor of two, it would be possible to give alert about 1 hour before coalescence. 
For an unequal-mass binary of (10 + 2 x 10^)Mq, the time the signal reaches p = 5 does 
not change much when we use adv. DPF (blue dotted). 

When only the ringdown signal of a galactic IMBH binary is detected by the 
ground-based detectors, DPF may help in determining the ringdown efficiency parameter 
erd = Ecv^/Mf [IHl HT] where -Eqw is the total energy emitted as radiation and Mf is 
the mass of the final BH. Since the ringdown amplitude is proportional to y/e^/Di |17] . 

and Dl degenerate. On the other hand, the direction and orientation of the binary 
can be determined fairly accurately due to the tremendous amount of SNR with the 
ground-based detectors. Then, from the inspiral and merger signals of DPF with SNR 
Pdpf, we can determine Dl with an error of AIuDl Pdpf- Using this information, we 
can determine from the ground-based interferometer data with an accuracy roughly 
the same as AIuDl. 

6. Probing the Mass of the Graviton 

In this section, we also consider the possible constraint on the mass of the graviton rUg 
from the GW observation with DPF. Originally, massive gravity theory was proposed 
by Fierz and Pauli [HI] where they simply added Lorentz-invariant mass terms to the 
Einstein-Hilbert action at the quadratic order. Nowadays, there are many kinds of 
massive gravity theories (see e.g. Refs. [95l [961 EH [HHl [99], llOOj ). Only recently a 
self-consistent massive gravity theory was proposed [99l llOOj which evades theoretical 
pathologies like Boulware-Deser ghosts under curved background [101] . In this paper, 
we do not stick any specific type of the massive gravity theories. All we assume is that 
the graviton has a finite mass nig. 

6.1. Previous Works for the Constraints on the Graviton Mass 

If a graviton has a finite mass, the gravitational potential becomes that of Yukawa 
type. This modifies the Kepler's third law and the current solar system experiment 
places the lower bound on the graviton Compton wavelength = h/{mgc) as Xg > 
2.8 X lO^^cm [51]. A weaker constraint has been obtained from the binary pulsar 
observations. Since there are additional gravitational degrees of freedom, the binary 
evolution due to gravitational radiation changes from that of general relativity and this 
can be tested from the orbital decay rate of the binary. By assuming the Fierz-Pauli- 
type theory, Finn and Sutton obtained the lower bound from PSR B1913+16 and PSR 
B1534+12 as > 1.6 x lO^^cm [52]. 

In the case of GW observations, the propagation speed of GW changes from the 
speed of light c and depends on its frequency / when the graviton is massive. This 
brings IPN correction in the phase of GW coming from a binary [102] . There are many 
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works that calculate the possible constraints on \g with future GW interferometers such 
as adv. LIGO, ET, LISA and DECIGO/BBO [1031 M, IM [1051 EB [271 HHl [T06l 
11071 1108] . We follow Keppel and Ajith |18] and use the hybrid waveform to estimate 
the lower bound on with DPF observations of IMBH binaries in our galaxy. 

6.2. Correction in the Gravitational Waveform Phase 

When the graviton has a finite mass m^, its phase velocity is modified from c as |102l [9T| 

(29) 

This leads to the correction to the phase of the gravitational waveform in the Fourier 
domain as |102[ US] 

^cfr(/) = ^(/)-/3r\ (30) 
where (3 is defined as 

Here, the distance D is different from the luminosity distance Di and it is given as |102] 

^-^r — . (32) 

Hq Jo (l + z')V^A/(l + ^')' + f^A 
For galactic IMBH binaries considered in this paper, the difference between Dl and D 
is very small and the factor [1 + z) can be approximated as 1 in the definition of (3 
(Eq. m)- 

6. 3. Fisher Analysis and Results 

We include /3 to the binary parameters and perform Fisher analysis to estimate 
how accurately we can determine the binary parameters (especially /3) with DPF 
observations. Then, we convert the upper bound on /3 to the lower bound on \g. First, 
we give a rough estimate on how strong we can constrain with DPF. It is not possible 
to detect the effect of finite graviton mass if the correction term in Eq. ( 1301) is 

smaller than SNR~^. This gives the constraint 



D 



A„ > 6.6 X 10"cm - — — — . (33) 

Next, we estimate the constraint numerically using Fisher analysis. In Fig. [SI we 
show the sky-averaged results for the lower bound on \g obtained for various BH masses. 
We assumed equal-mass binaries with x = 0.2 at Dl = 5kpc. The results using DPF are 
shown in the (red) thick solid curve. For a (10^ + 10^)Mq binary, the constraint becomes 
Xg > Qx lO^^cm. This is weaker than our rough estimate in Eq. (|33l) and this is due to the 
degeneracies between /3 and other binary parameters. The results using adv. DPF are 
shown in the (green) thin solid curve. We can understand this behaviour by combining 
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Figure 8. The lower bounds on Ag using DPF with (red thick sohd) and without 
(magenta thick dotted) taking x into parameters, and using adv. DPF with (green thin 
sohd) and without (blue thin dotted) taking x into parameters. We assumed equal- 
mass binaries with x = 0-2 at D^, = 5kpc. The (black) thin dotted horizontal line at 
Xg = 2.8 X lO^^cm represents the bound obtained from solar system experiment [5T| 
while the (black) thick dotted-dashed horizontal line at Ag = 1.6 x lO^^cm shows the 
one obtained from binary pulsar test |52) . 



Fig. [Hand the fact that the constraint on is proportional to SNR~^/^ (see Eq. f p3|) ). 
In this paper, we only focus on binaries with aligned spins. When we consider more 
general cases, there will be precessions which solve the degeneracies between spin and 
other parameters. Just to give some idea for the results in these situations, we calculated 
the constraints on without taking x as a variable parameter. (See Ref. |105] for the 
discussion that the constraint on A^ including precession is almost the same as the 
one without taking spins as variable parameters.) This corresponds to the situation 
where the degeneracies between A^ and x are completely solved. The results using 
DPF and adv. DPF are shown in the (magenta) thick and the (blue) thin dotted curve, 
respectively. The dotted-dashed horizontal line at A^ = 2.8 x lO^^cm corresponds to the 
(static) lower bound obtained from the solar system experiment [51]. Although DPF 
constraint is slightly weaker than this, it is still meaningful since the DPF measures 
the deviation in the propagation speed of GWs while the solar system experiment 
measures the deviation in the gravitational constant (or in the Kepler's third law). 
Our results show that DPF can put about 2 orders of magnitude stronger (dynamical) 
constraint than the one obtained in the weak-field test of the binary pulsar [52] which 
is shown as the dotted-dashed horizontal line at A^ = 1.6 x lO^^cm. Furthermore, Finn 
and Sutton [52] assumed Fierz-Pauli-type theory while the constraint obtained here is 
independent of the specific massive gravity theory. Since D^^ in the amplitude cancels 
with D in f3 when calculating the Fisher matrix (see Eq. (l33l) ). the results estimated 
here are almost independent of Dl- The constraint is stronger for larger mass binaries, 
hence the constraint becomes stronger if we can also reduce the acceleration noises. 
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Figure 9. Histograms sliowing the lower bounds on Ag obtained from Monte Carlo 
simulations with lO^M© equal-mass IMBH binaries in Centauri with random 
orientation for the orbital angular momenta. We show here the results using DPF. 
The (black) thick and (red) thin histogram each represents the one with and without 
taking % as a variable parameter, respectively. 

Next, we performed Monte Carlo simulations for IO^Mq equal-mass binaries of 
X = 0.2 in uj Centauri. The results using DPF are shown in Fig. |9l The (black) 
thick histogram and the (red) thin one each shows the probability distribution for the 
upper bound on \g with and without taking x as a variable parameter, respectively. 
This shows that the constraint becomes slightly weaker compared to the sky-averaged 
analysis. However, again, this is much stronger than the one from binary pulsar tests. 
When we use adv. DPF, these probability distributions shift to larger \g by roughly 
v^ = 2.6. 

7. Event Rate Estimations 

7.1. Mergers of Equal-Mass IMBH Binaries Formed in Galactic Massive Young 
Glusters 

Currently, more than 10 galactic massive young clusters (GMYCs) have been 
discovered [109j . Giirkan et al. |13] performed numerical simulations and found that 
GMYCs may contain two IMBHs at their centers with BH masses IO^Mq, which are 
likely to form binaries. After IMBH binary formation, it shrinks due to the dynamical 
friction with the cluster stars. The timescale of this process is typically < 10 Myr which 
is independent of the local average stellar- mass [l3l IllOj . Then, binary shrinks via 
dynamical encounters with cluster stars, with the timescale of < 1 Gyr |110] . Finally, 2 
IMBHs merge due to gravitational radiation within 1 Myr [43]. 

For simplicity, we assume that IMBH binaries are all situated at the galactic center. 
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Table 2. The distances and velocity dispersions of galactic globular clusters. Possible 
masses of IMBHs, if they exit, are obtained from M — a relation [112) . 



NGC 


distance 


vel. disp. a 


BH mass 


No. 


(kpc) [63] 


(km/s) [mj 


(Mo) 


104 


4.5 


10.0 


794.7 


362 


8.5 


6.2 


116.3 


1851 


12.1 


11.3 


1299 


1904 


12.9 


3.9 


18.04 


5272 


10.4 


4.8 


41.57 


5286 


11.0 


8.6 


433.4 


5694 


34.7 


6.1 


108.9 


5824 


32.0 


11.1 


1209 


5904 


7.5 


6.5 


140.6 


5946 


10.6 


4.0 


19.97 


6093 


10.0 


14.5 


3539 


6266 


6.9 


15.4 


4508 


6284 


15.3 


6.8 


168.6 


6293 


8.8 


8.2 


357.9 


6325 


8.0 


6.4 


132.4 


6342 


8.6 


5.2 


57.35 


6441 


11.7 


19.5 


11645 


6522 


7.8 


7.3 


224.3 


6558 


7.4 


3.5 


11.68 


6681 


9.0 


10.0 


794.7 


7099 


8.0 


5.8 


88.96 



From Fig. |3l we see that DPF is not sensitive enough to detect GW signals from IMBH 
binaries in GMYCs. On the other hand, when we use adv. DPF, we will be able to see 
all of the IMBH binaries in GMYCs. Following Fregeau et al. |110] . we assume that 
the number of star clusters massive enough to form IMBH binaries equals to the one 
of globular clusters, and 10% of them actually produce IMBH binaries. In our galaxy, 
there are about 150 globular cluster [63], hence 15 IMBH binaries are expected to be 
within the reach of DPF. Since only one IMBH binary is formed over its lifetime for 
each cluster jllO] , the detection rate of the IMBH binaries in GMYCs with adv. DPF 
can be roughly estimated as z^gmyc ~ 15/(13. 8Gyr) = 1.1 x 10~^yr~^. 

7.2. Equal-Mass IMBH Mergers in Globular Glusters 

IMBH binaries at the center of galactic globular clusters are also potential GW source 
candidates for DPF. Based on Table 6 of Ref. jlllj . we select 21 globular clusters out 
of 150 globular clusters that have been found. Since the velocity dispersion (denoted as 
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Figure 10. Same as Fig. |3l but here we also plot possible IMBH binaries shown in 
Table [2] as " x " with masses assumed from M — a relation. 



(Tp(core)) are listed in the table, we apply the M — a formula obtained by Tremaine et 
al. [m] as 



to estimate the mass of the assumed IMBH in each globular cluster. We assume that 
this BH is composed of an equal-mass IMBH binary with each BH having mass M/2. 
The results are listed in Table [2], together with the distance to each globular cluster 
which is obtained from the catalogue of galactic globular clusters made by Harris (see 
the URL link shown in Ref. [63]). We plot IMBH binary of each globular cluster in 
Fig. [To] and count how many of them have SNRs larger than 5 so that they can be 
detected by DPF. The meaning of this figure is same as Fig. [31 We see that about 2 out 
of 26 globular clusters (5 shown in Table [U + 21 shown in Table W) may contain IMBH 
binaries that might be detected by DPF on average. Then, the detection rate of the 
IMBH binaries in globular clusters vqc is given as 

150 1 in 1 

i/GC ~ 2 X X — = 8.4 X 10"^°yr-\ 35 

26 13.8 Gyr ^ ^ 

Here, we divide the number of globular clusters by the age of the universe for the same 

reason discussed in the previous subsection. When we use adv. DPF, we will be able to 

see 13 out of 26 IMBH binaries in the galactic globular clusters, making the event rate 

larger than above by a factor 13/2 = 6.5. 

7.3. Intermediate-Mass Ratio Inspirals (IMRIs) with Adv. DPF 

Many IOMq BHs are expected to form after supernova explosions and they sink to 
the cluster cores due to the mass segregation. Therefore at the centers of the clusters, 
these BHs may have number densities comparable to the main sequence stars, typically 
n = lO^pc"^. In this subsection, we consider IMBH binaries with intermediate-mass 
ratio inspirals (IMRIs). In Fig. [TH we show observable ranges for DPF (red thick solid) 
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Figure 11. Observable ranges against total mass M with m2 = lOM©. The meaning 
of each curve and plot is same as Fig. [TOl 



and for adv. DPF (green thin solid) with m2 = IOMq. Meanings of curves and plots are 
same as in Fig. [TOl We can see that DPF cannot observe any IMRIs with m2 = IOMq 
but adv. DPF is sensitive enough to possible IMRIs in some of GCs and may be able to 
detect an IMRI signal from the galactic center. 



7.3.1. Maximum Merger Rate Per Cluster First, following Ref [50], we will estimate 
the upper bound for the merger rate per cluster, which is set from (i) the supply limit 
of smaller BHs, and (ii) the mass of the larger BH. For case (i), the rate must not 
exceed z/guppiy = A^Bn/^age where N^h is the current number of smaller BHs and Tage 
is the current age of the cluster. Otherwise, almost all of the smaller BHs have been 
consumed by now and A^bh would have been reduced considerably. 

For globular clusters with mi = IO^Mq, let us assume that the number of stellar- 
mass BHs with m2 = IOMq is 10'^ (this value is used for NSs with central BH mass 
nil = IO^Mq in Ref. [80]). Then, the rate for the supply limit becomes 

For galactic massive young clusters (GMYCs) with ages less than 3 x lO'^'yr and 
with masses larger than IO^Mq, the number of main sequence stars are at least 
10^Mq/0.4Mq = 2.5 X 10^. For a multimass King model, the scale height of the smaller 
BH m2 is smaller compared to the main sequence stars mms by a factor (m2/mms)^/^ [65] . 
This means that when the number densities of the main sequence stars and the smaller 
BHs are the same at the cluster core, the number of smaller BHs is smaller by a factor 
{m2/m^sY^'^ = 125 (we here assumed that m2 = IOMq and = 0.4Mq) compared 
to the main sequence stars. Therefore, the number of smaller BHs is roughly 2 x 10^ 
and the supply limited rate becomes 

GMYC _ 2 X 10^ _ 7 . n-Vr-l C^r) 
'^supply - 3x lOV ~ ■ 
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Figure 12. The SNR distribution for lO^M© equal -mass IMBH binaries in young 
clusters at the galactic center. We randomly distribute the orbital angular momentum 
orientations. The (red) thin and the (black) thick histogram each represents the one 
for the static and the moving detector, respectively. The vertical dashed line at p = 5 
denotes the SNR threshold. 



For case (ii), the merger rate cannot exceed t'mass = ("^i/''^2)^ago since otherwise, 
the mass of the larger BH has been increased considerably. For globular clusters with 
mi = 1O^M0, this limit is set as 



GC 



WMq 1 

IOMq 13.8Gyr 



while GMYCs give 

,GMYC 







1 



7.2 X 10"V~\ 



3.3 X 10~V~^- 



IOMq 3 X lOV 
In total, the upper bound of the merger rate is set as 

^\im\t min [l^supplyi ^mass] ^mass' 



(38) 

(39) 
(40) 



7.3.2. Merger Rate Next, we estimate the merger rate of IMRIs in GCs and GMYCs, 
with the limit discussed in the previous subsection. For GCs, we use t'l^mit — ^mass ^ 
our merger rate per cluster lffl . From Fig. [TTl we see that IMRIs in about 3 out of 26 
globular clusters can be seen by adv. DPF. This leads to the estimate that it can see 
IMRIs from about 3 x 150/26 = 17.3 globular clusters in total, hence the detection rate 
of globular clusters in total becomes 

^GCtot < 1.2 X 10- V~^- (41) 

tfHopman [TT^ estimated the merger rate of IOMq BH into 3 x IO^Mq IMBH with Monte Carlo 
simulations, obtaining v k. T x 10~^ yr"^. If we assume that the order does not differ much for smaller 
IMBH mass cases, this exceeds the limit in Eq. (|40p . 
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For GMYCs, from Fig. [TTl we see that adv. DPF is not sensitive, on average, to an 
IMRI signal of (10 + 10'^)Mq at galactic center. However, for some binaries with near 
optimal orientations of the angular momentum X, SNRs go beyond p = 5 with DPF. 
Let us consider what fraction of these binaries can be seen by DPF (for simplicity, we 
here fix the BH mass as lO^Mo). From Eq. (ITTl) . the angle dependence of SNR squared 
is given as 

p2 oc oc sin^ ^{cos^ 2ip{l + cos^ if + 4 sin^ 2^ cos^ i} = X. (42) 

For a (10 + 10^)Mq IMRI at Dl = 8.5kpc with p = 5, X becomes X = Xq = 3.13. 
Then, the probability that X being greater than Xq is 

P{X>Xo) = (^0 dfie j dfi,j^ dip 

e[(l - plficos" 27/>(l + /i2)2 + 4 sin2 2^p2| _ 
= 0.027, (43) 

where pe = cos 9, ji-i = cosz and 6[- ■ ■] represents the step function. 

In order to confirm this analytic estimate, we performed Monte Carlo simulations. 
We first consider static detectors with (fiEit) = v^eo and (pB(t) = V^do- Then, we randomly 
distribute 10"^ sets of cosi9l in the range [-1,1] and 0l, V^eo and ipDo in the range [0,27r]. 
We set the galactic latitude and longitude of the source as (6, /) = (0, 0) (galactic center). 
We calculate SNR of each set and count the fraction of binaries P{X > Xq) that have 
SNR greater than 5. The result of probability distribution function is shown in Fig. [12] 
with (red) thin histogram. The vertical dashed line corresponds to p = 5. In this case, 
we found P{X > Xq) = 0.028 which shows good agreement with Eq. ( l43i) of our analytic 
estimate. Next, we take the detector motion into account. We take (pE(t) = <^Et + V^eo 
and </?d(^) = + v^DO; and performed the same Monte Carlo simulations. The (black) 
thick histogram in Fig. [12] is the one obtained with moving detectors. This time, the 
fraction becomes P{X > Xq) = 0.017. When we take the detector motion into account, 
SNR approaches to the sky-averaged value and the probability distribution becomes 
sharper. This is why P{X > Xq) appeared smaller than the one for the static detector 
case. If we manage to improve the sensitivity of adv. DPF twice, SNRs would be doubled 
so that the probability of IMRIs with SNR greater than 5 becomes roughly 50%. This 
means that the detection rate would be improved by about 30 times. 

According to Ref. |114] . the merger rate due to three-body interaction per 
cluster becomes i/G;^yc _ 2 x 10~*yr~^. (For relatively large IMBH masses, the 
two-body interaction may give comparable contribution compared to the three-body 
interactions [80].) Notice that this value does not exceed the limit we found in Sec. 17.3.11 
Gvaramadze et al. [109] proposed that there should be 70-100 massive young clusters 
in the galactic disk in total. They also suggested that about 50 young massive clusters 
exist in the galactic center. Following Ref. [114] , we assume that 75% of GMYCs contain 
IMBHs at their centers, meaning that about 50-75 GMYCs may contain IMBHs. We 
make a conservative assumption that 50 GMYCs may contain IMBHs. In this case, the 
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total detection rate becomes 



,GMYC,tot 



0.017 X 50 X 2 X lO^V 



1.7 X lO^V 



-1 



(44) 



8. Conclusions and Discussions 

DPF has an ability to detect IMBH binaries within our galaxy. Although there is plenty 
of evidence for the existences of stellar-mass BHs and SMBHs, there is no direct evidence 
of an IMBH and its existence is still controversial. IMBHs may exist at the centers of 
the galactic globular clusters and the galactic massive young clusters (GMYCs). For 
the former cases, globular clusters such as u Centauri |36], M15 [35] and NGC 6752 [IQ] 
may harbor IMBHs with masses larger than IO^Mq. For the latter, it is suggested that 
there may exist about 100 GMYCs which contain IMBHs |109j . A numerical simulation 
shows that these GMYCs may contain IMBH binaries with WMq BHs [13]. 

If a IO^Mq equal-mass IMBH binary exists in cu Centauri, DPF may perform unique 
observation since this mass range is too large for the ground-based detectors. On the 
other hand, if a lO^M© equal-mass IMBH binary exists in NGC 6752, DPF can see 
the late inspiral part while LCGT may detect merger and ringdown signals. Ringdown 
signal may also be detected by adv. LIGO and adv. VIRGO. 

We performed both pattern-averaged analysis and Monte Carlo simulations to 
estimate how accurately we can determine the binary parameters when the signal is 
detected. For the latter simulations, we take the motion of DPF into account. If the 
one from uj Centauri is observed, it seems that the total mass can be determined with 
a 2-3% accuracy, while the symmetric mass ratio and the effective spin parameter can 
be determined with 5-6% accuracies. 

We also discussed the possible contributions of DPF to the ground-based GW 
searches. First of all, DPF may give alert to the ground-based detectors about 10 
mins before coalescence for 10"^ Mq equal-mass IMBH binaries in NGC 6752. Therefore, 
the ground-based ones can get ready for their ringdown detections. Also, DPF data 
including inspiral and merger signals of IMBH binaries may help in distinguishing the 
ones from the noises when only ringdown signals have been detected on the ground. 
Furthermore, by combining DPF and the ground-based data, we may be able to 
determine the ringdown efficiency which cannot be obtained from the ringdown 
data alone due to the degeneracy between erd and the distance to the BH. 

Also, it may be possible to constrain the graviton Compton wavelength A^, with 
slightly weaker bound compared to the solar system result in the weak- field regime |51j . 
DPF constraint, independent of the massive gravity theory, is almost 2 orders of 
magnitude stronger than the binary pulsar tests [52], which has been obtained by 
assuming Fierz-Pauli-type theory [M] . 

Unfortunately, the detection rate of DPF is not so high. For equal-mass IMBH 
binaries, it is roughly estimated as 10~^ yr"^- However, when we use adv. DPF, it has 
an ability to observe IMBH binaries with m2 = IOMq which is a typical mass for the 
remnant BH of supernova and due to the mass segregation, it is expected that there 
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are plenty of IOMq BHs at the centers of clusters. Therefore, the detection rate would 
be larger for these binaries. For example, the one of (10 + 10'^)Mq in massive young 
clusters at the galactic center is estimated as 10~* yr~^-10~'' yr~^. Given these values, 
the prospects of detecting IMBH binaries or IMRIs with IMBHs using DPF or adv. DPF 
are not very optimistic. 

In this paper, we considered the possibility of reducing the laser frequency noise 
down to the shot noise. It is also interesting if we can reduce the acceleration noises as 
well. This has more impact on larger mass binaries (like the one in u Centauri that we 
assumed in this paper) and, for instance, the constraint on becomes stronger with the 
scaling proportional to SNR^/^. Also we will be able to give earlier alert to the ground- 
based detectors if we can improve the sensitivity on lower frequency part. Although the 
expected detection rate is rather small, the possibility is not zero and since we would 
be able to obtain interesting science from these binaries (that we discussed in the main 
part of this paper), we highly recommend to reduce the noise and accomplish at least 
the sensitivity of adv. DPF. 
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Appendix A. Relationship between the Celestial Coordinate and the 
Galactic Coordinate 

We first introduce the right ascension a and the declination 5 shown in the panel (a) of 
Fig. lAll The direction of the source (^s, 0s) can be expressed as {6a,, 0s) = (f ~ o;)- 
Next, we introduce the galactic coordinate with the galactic latitude h and longitude / 
shown in the panel (b) of Fig. lAll (6, /) is related to (a, 5) as [115] 



where (aGP,^Gp) = (192.86°, 27.13°) represents the North Galactic Pole (NGP) and 
lev = 123.93° is the longitude of the North Celestial Pole (NGP). From these equations, 



sin 6 = sin 6qp sin b + cos 5gp cos b cos(/cp — ; 

cos5sin(a — acp) = cos6sin(/cp — /), 

cos 6 cos(a — ctGp) = cos 6gp sin b — sin 6gp cos b cos(/cp — 



(A.l) 
(A.2) 
(A.3) 



(^S) 0s) can be written as 



cos 9s = sin 6gp sin b + cos 6gp cos b cos(/cp — > 



(A.4) 
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Figure Al. (a) The celestial coordinate {a;,y, z} with the right ascension a and the 
declination 5. x and y-axes exist on the celestial equator with the former pointing 
towards the vernal equinox, z points towards the North Celestial Pole (NCP). (b) The 
relation between the celestial plane and the galactic plane with the galactic latitude 
h and longitude I. b and I are measured from the galactic equator and the galactic 
center, respectively. The galactic equator is tilted 62.5° from the celestial one. NCP 
stands for the North Galactic Pole. 

sin^s = a/1 — {sin^GP sin 6 + cos^gp cos6cos(/cp — /)}^, (A. 5) 

cos 0s = [ — cos h sin(/cp — I) sin ogp 

+ {cos^GP sin 6 — sin^GP cos6cos(/cp — 0} cosc^gp]/ cos 5, (A. 6) 
sin 0s = [ cos h sin(/cp — I) cos ogp 

+ {cos^GP sin 6 — sin^GP cos6cos(/cp — 0} sinc^Gp]/ cos 5. (A. 7) 
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